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ABSTRACT 

V2491 Cygni (Nova Cygni 2008 No. 2) was detected as a transient supersoft X-ray source with the Swift XRT 
as early as 40 days after the outburst, suggesting a very massive white dwarf (WD) close to the Chandrasekhar 
limit. We present a unified model of near infrared, optical, and X-ray light curves for V2491 Cyg, and have 
estimated, from our best-fit model, the WD mass to be 1.3 ± 0.02 M Q with an assumed chemical composition 
of the envelope, X = 0.20, Y = 0.48, X CNO = 0.20, X Ne = 0.10, and Z = 0.02 by mass weight. We strongly 
recommend detailed composition analysis of the ejecta because some enrichment of the WD matter suggests 
that the WD mass does not increase like in RS Oph, which is a candidate of Type la supernova progenitors. 
V2491 Cyg shows a peculiar secondary maximum in the optical light curve as well as V1493 Aql and V2362 
Cyg. Introducing magnetic activity as an adding energy source to nuclear burning, we propose a physical 
mechanism of the secondary maxima. 

Subject headings: novae, cataclysmic variables — stars: individual (V1493 Aql, V2362 Cyg, V2491 Cyg) — 
stars: mass loss — X-rays: binaries 



1. INTRODUCTION 

Classical novae show a wide variety of timescales and 
shape s in the optical light curves (e.g., Pavne-Garjoschkin 
fl957h . Among various shapes of nova light curves, V1493 
Aql (Nova Aquilae 1999 No.l) shows an impressive sec- 
ondary maximum about 50 days after the outburst (e.g., 
iBonifacio et al. 1 120001: IVenturini etaTl l2004h . although the 
physical mechanism of the secondary maximum is not under- 
stood yet. Recent two novae, V2362 Cyg (Nova Cygni 2006) 
and V2491 Cyg (Nova Cygni 2008 No. 2), also show a similar 
type of single secondary maximum, at about 250 and 15 days 
after the outburst, respectively. These three novae form a wide 
variety set of timescales, i.e., about 15, 50, and 250 days at the 
secondary maximum and of secondary p eak heights, i.e., 1.1, 
2.8, a nd 3.6 mag, respectively, (see, e.g jKim eswenge r et al. I 
2008), which provide us a new clue to the mechanism of the 
secondary maxima. 

In this Letter, we propose a strong magnetic activity as the 
mechanism of the secondary maxima observed in V2491 Cyg, 
V1493 Aql, and V2362 Cyg, using the white dwarf (WD) pa- 
rameters obtained from light curve fittings based on an op- 
tically thick wind model of nova outbursts dKato & Hachisu I 
1994). In §2, we briefly describe our numerical method and 
light curve fitting of V2491 Cyg. In §3, we propose an idea of 
strong magnetic activity in the WD envelope and estimate the 
timescales of the secondary maxima. V1493 Aql and V2362 
Cyg show very different timescales of the secondary maxi- 
mum, both of which are also explained by the same mecha- 
nism in §4. Conclusions follow in §5. 



2. MODELING OF NOVA OUTBURSTS 
2.1. Optically thick wind model 
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After a thermonuclear runaway sets in on a mass-accreting 
WD, its photospheric radius expands greatly to R p h > 100 Rq 
and the WD envelope settles in a steady-state. We have fol- 
lowed evolutions of novae by connecting steady state solu- 
tions along the decreasing envelope mass sequence. We solve 
a set of equations, that is, the continuity, equation of motion, 
radiative diffusion, and conservation of energy, from the bot- 
tom of the hydrogen-rich envelope through the photosphere 
assuming spherical symmetry. Winds are accelerated deep in- 
side the photosphere so that they are called "optically thick 
winds." As one of the boundary conditions for our numeral 
code, we assume that photons are emitted at the photosphere 
as a blackbody with the photospheric temperature of 7J,h. X- 
ray flux is estimated directly from the blackbody emission, but 
infrared and optical fluxes are calculated from free-free emis- 
sion by using the physical values of our wind solutions. We 
neglect the effect of ash helium layer, which may be piled up 
beneath the hydrogen bu rning zone, for all nov a calculations 
except the RS Oph case (Hac hisu et al. Il2007h . Our method 
and various physical proper ties of these wind solutio n s have 
already been published (e.gjHachisu & Kato 2001a b j2004 , 
20061 [2007: Ha chisu et all 1 9961 11999a b. 2000, 2003,|2007, 
2008; K atoll 1 983L 1 1 997l 1 1 999fc iKato & Hachisu 111994 . 

The light curves of our optically thick wind model are pa- 
rameterized by the WD mass (MwdX chemical composition 
of the envelope (X, Y, Xqno, X^ e , and Z), and the envelope 
mass (AM env .o) at the outbur st (day 0). Detail s of our ligh t 
cur ve fittings are described inlHachisu & Ka to (2006, [2007h 
and lHachisu et al. I (I2007ll2008l) ~ 

2.2. Light curve fitting (V2491 Cyg) 

V2491 Cyg was discovered by Nishiy ama and Kabashim a 
at mag 7.7 on 2008 April 10.728 UT (Nak ano et al~ll2008h . 
The nova was not detected on April 8.831 UT (limiting mag 
14). The exact outburst day is unknown, so we assume here 
that ?ob = 2454566.0 (April 9.5 UT) is the outburst day (day 
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V2491 Cyg (Nova Cyg 2008#2) 




log[days after outburst] 

FIG. 1 . — Model light curves of V2491 Cyg for three white dwarf masses, 
1.32 Mq (green line), 1.3 M Q (thick bl ue line), and 1.27 Mp, ( thin black 
line), together with the observations. See Hachisu & Kato 12006) for more 
details of light curve fitting. Two arrows indicate epochs when the wind 
stops (day 40) and when the hydrogen shell-burning ends (day 50) in the 
1.3 Mq WD model. Dashed lines: Supersoft X-rays are probably not de- 
tect ed during the wind phase because of self-absorption by wind itself (see, 
e.g. [Hachisu & Kato 2003b). Large open circles at the right end of each op- 
tical light curve denote the epoch when the optically thick wind stops. We 
obtain the best-fit model for the envelope chemical composition of X = 0.20, 
Y = 0.48, X CN0 = 0.20, X Nc =0.10, and Z = 0.02. Large open tri angles: Ob- 
serva tional X-ray (0.2—0.6 keV) count rates obtained with Swift I P age et al. I 
2009). Optical and near IR observational data of / (open diamonds), R (cir- 
cles with a plus), V (filled squares), y (filled circles), and visual (small open 
circles) are taken from AAVSO (American Association of Variable Star Ob- 
servers) and VSOLJ (Variable Star Observers League in Japan). The F\ oc f~ 3 
law (magenta) is added for the nebular phase, i.e., after the wind stops, where 
1 is the time after the outburst. 

0). The orbital period of P OT b = 0.0958 days was derived by 
iBaklanov et al. I d2008l) from the modulations with an ampli- 
tude of 0.03 - 0.05 mag. 

The rise or decay time of supersoft X-r ay flux is an im- 
portant indicator of the WD ma ss (e.g., lHachisu & Kato I 
120061 120071; lHachisu et al71l2008l) The best fit model with 
the Swift observation dPage et al. 1120091) is the WD mass of 
1.3 ±0.02 Mq as shown in Figure Q] Here we adopt the 
chemical composition of X = 0.20, Y = 0.48, Xcno = 0.20, 
X Ne = 0.10, and Z = 0.02 to reproduce a short (10 days) 
supersoft X-ray duration ( see lHachisu & Kato I I2006L 12007 ; 
lHachisu et al. 1120071 120081 for dependence on chemical com- 
position). 

Figure Q] also shows optical and near IR light curves of the 
/, R, V, y, and visual magnitudes. Our theoretical light curves 
of free-free emission reasonably fit with the observation until 
day ~ 50 except the secondary peak. After that, the observa- 
tional data deviate probably due to the contribution of nebular 
emission lines. 

Thus we conclude that the WD of V2491 Cyg is as mas- 
sive as 1 .3 Mq from our light curve fittings with the supersoft 
X-ray, optical, and near IR data. This WD mass of V2491 
Cyg is a bit less mas sive but comparable to that of the re- 
current nova RS Oph (lHachisu et al. Il2007l) . the supersoft X- 
ray duration of which lasted 60 days startin g from day 30, 
much longer than the 10 days in V2491 Cyg. Hachis uet al. I 
d2007l) explained this long duration of supersoft X-ray phase 
of RS Oph by incorporating heat flux from a hot helium layer 
underneath the hydrogen burning zone, which develops in 
mass-increasing WDs. In V2491 Cyg, however, the supersoft 
X-ray phase lasts only 10 days, indicating that no thick he- 
lium layer develops beneath the hydrogen burning zone and, 
therefore, the WD m a ss is not increasing (see discussion in 
lHachisu et al. 112007b . iTomov et al. I (120081) suggested a pos- 
sibility that V249 1 Cyg is a recurrent nova, for which we ex- 
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FIG. 2. — Five energy densities of the WD envelope are plotted against the 
radius, i.e., thermal energy (e^: red), thermal gas energy (e t h.gas : magenta), 
rotational kinetic energy (e ro t: blue), magnetic energy (emag: black), and 
wind kinetic energy (£ w i„d: green). Five sequential stages are plotted for e m t 
but only two stages [the first (thick solid) and the last (thin solid)] for the other 
energy densities. Open circles indicate the photosphere while open squares 
correspond to the critical point of each wind solution. Various parameters are 
summarized in the figure. 

pect that the WD mass is increasing ( lHachisu & Kato 12 00 lb). 
One of the strongest clues to this question is the chemical 
composition of the ejecta. We strongly recommend compo- 
sition analysis of V2491 Cyg to clarify whether or not the 
ejecta are conta minated by WD ma tter. Suzaku X-ray spec- 
tra analyzed by Tak ei et al71 ([2009 ) suggested oxygen/neon- 
rich ejecta. If so, the WD mass is decreasing and, as a result, 
V2491 Cyg is not a candidate of Type la supernova progeni- 
tors. 

3. MAGNETIC ACTIVITY (V2491 CYG) 

V2491 Cyg shows a single secondary maximum in the op- 
tical and near IR light curves, which we cannot reproduce by 
our evolution model. Here we propose magnetic activities as 
a new energy source of the secondary maximum. 

The prenova X-ray was de tected in V2491 Cyg with Swift 
(llbarra & Kuulkers 1 120081) . Ilbarra et al. I (12008b suggested 
that the prenova X-ray spectra are more like those seen 
from magnetic ra t her th an non-magnetic cataclysmic vari- 
ables. iTakei et al71 (12009) reported the first detection of super- 
hard X-rays (15-60 keV) with the Suzaku HXD at day 10, the 
spectrum of which has a power law distribution, suggesting a 
non-thermal origin, and unlikely to arise from thermal emis- 
sion from internal shocks. This superhard component was not 
seen in their second Suzaku observation at day 30. These two 
observations point toward strong magnetic activities on the 
WD surface. 

Therefore, our idea on the secondary maximum is based on 
additional energy release associated with rotating magnetic 
field. We assume that V2491 Cyg is a polar system with mag- 
netic field as s trong as Z?o ~ 10 7 G on the WD surface (e.g., 
IWarner 111 9951) . Before the nova outburst, the WD magnetic 
field rotates synchronously with the WD spin as well as the 
binary orbital motion. After the outburst, the nova envelope 
expands largely and rotates differentially due to local angular 
momentum conservation. Then, the magnetic field no more 
rotates synchronously with the WD spin because the magnetic 
tension is not strong enough to keep the whole envelope rotat- 
ing with the WD spin. We expect that the differential rotation 
amplifies magnetic field which drives strong magnetic activi- 
ties. As the nova outburst proceeds, the density of the WD en- 
velope is gradually decreasing due mainly to wind mass loss 
and then the magnetic field eventually recovers synchronous 
rotation. The strong magnetic activities end at this stage, cor- 
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TABLE 1 

Epoch of Secondary maximum in our nova models 



object 


max" 
(day) 


(day) 


M 2 b 
(Mq) 


a 

(«©) 


M W D 

(M ) 


£rot ~ £mag 

(day) 


V2491 Cyg 


15 


0.0958 


0.18 


1.0 


1.32 


13 












1.3 


15 












1.27 


18 


VI 493 Aql 


50 


0.156 


0.34 


1.4 


1.2 


40 












1.15 


49 












1.1 


62 


V2362 Cyg 


250 


0.207 


0.48 


1.6 


0.75 


200 



0.7 240 

0.65 330 

a epoch of secondary maximum^ calculated from equation (5} 

responding to the end of a secondary maximum. Since the 
additional energy source disappears, the light curve goes back 
to a "normal" one as shown in Figure[T] 

To confirm this idea we estimate the total thermal energy 
of gas plus radiation (e t h = £th,gas + £th,rad), rotational kinetic 
energy e rot , wind kinetic energy e w i n d, and magnetic energy 
£ma g > which are calculated from 

3 kT 

Eth = Eth.gas + Efh.rad = T p+ClT 4 , (1) 

2 limn 

1 2 
£rot= 2/°( r ^spin) , (2) 

_ 1 2 

£\vind — ^^windi (3) 

«--£Gs;f- (4) 

The temperature T, density p, and wind velocity v w j n d are 
taken from our wind solutions of the best-fit model of V249 1 
Cyg (1 .3 Mq WD). Here we assume the dipo le magnetic f ield 
of B = 3 x 10 7 G at the WD surface (e.g.. lWarnerlll995L for 
VI 500 Cyg) and that the WD spin period is the same as the 
orbital period, 27r/S7 sp i n = f sp i n = P orD . Figure [2] shows the 
energy densities for five sequential stages during the nova 
outburst. The first model (thick solid) has the largest pho- 
tospheric radius and corresponds to a stage at/near the op- 
tical maximum. We see that at/near the optical maximum, 
e mag < e rot at r > 2 x 10 9 cm. This indicates that the magnetic 
field is differentially rotating outside of r ~ 2 x 10 9 cm. Then, 
the rotation energy density decreases with time and eventually 
£mag becomes comparable to e rot at the stage of the smallest 
photospheric radius. After that, the magnetic field probably 
gains synchronous rotation with the WD spin. We expect that 
the magnetic activities have a peak at e mag rj e rot . This con- 
dition is satisfied at day 15 for our best-fit model, being very 
consistent with the time of the secondary maximum. 

Next we estimate the epoch when the companion emerges 
from the nova envelope. If the mass of the donor star is esti- 
mated from Warner's (1995) empirical formula, i.e., 

^ ~ 0.065 f-^V /4 , for 1.3 < < 9 (5) 
M© \ hours J hours 

we have = 0.18 Mq, which corresponds to the separa- 
tion of a = 1.0 Rq, and the effective Roche lobe radius of 
the primary component (WD) of = 0.55 R Q . When the 
photospheric radius of the nova envelope shrinks to near the 
orbit (the separation a), the condition of e rot < e mag is sat- 
isfied as shown in Figure [2] This indicates that the epoch 



V1493 Aql (Nova Aql 1 999#1 ) 
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FIG. 3.— Same as Fig.Q] but for V 1493 Aql (Nova Aql 1999 No.l). Our 
best-fit model is Mwd = 1.15 Mq for the envelope chemical composition of 
X = 0.55, ¥ = 0.30, Xctm = 0.10, X Ne = 0.03, and Z = 0.02. Observational 
data are taken from AAVSO and VSOLJ. We further add the data of IAU 
Circ. 7223, 7225, 7228, 7232, 7258, 7273, 7313. 
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FIG. 4.— Same as Fig. \T\ but for V2362 Cyg (Nova Cyg 2006). Our 
best-fit model is A/wd = 0.7 Mq for the envelope chemical composition of 
X = 0.45, ¥ = 0.18, X CN0 = 0.35, and Z = 0.02. Observational data are taken 
from AAVSO and VSOLJ. 

of the secondary maximum is shortly after the companion 
emerges from the nova envelope. We here implicitly assume 
that strong magnetic field connects the WD and the compan- 
ion (like in polar systems). The mechanism of activity we sup- 
pose is magnetic reconnection. Strong magnetic reconnection 
occurs between the WD and the companion. When the WD 
photosphere is larger than the companion's orbit, magnetic re- 
connection may occur deep inside the photosphere but the gas 
pressure (or gas thermal energy) at the reconnection region is 
much larger than the magnetic pressure (magnetic energy), so 
the gas is not easily accelerated by magnetic force. On the 
other hand, when the gas pressure becomes smaller than the 
magnetic pressure, i.e., when the photosphere shrinks to the 
orbit (see Figf2]l, then the gas is easily accelerated by mag- 
netic force and the envelope gas is massively ejected. Thus 
this process increases the mass-loss rate around/near the sec- 
ondary maximum. 

4. V1493 AQL AND V2362 CYG 

V1493 Aql was discovered by Tago at mag 8.8 on July 
13.558 UT. The nova was not detected on his films of July 
5 and 9 (limiting mag 11 and 10.5, respectively). There- 
fore, we assume here that the outburst day is ?ob = 2451372.0 
(July 12.5 UT). Figure[3]shows optical and near infrared light 
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curves of V1493 Aql. Our theoretical light curves of free-free 
emission reasonably fit with the observation until day ^100 
except the secondary maximum. The best-fit model is the 
WD mass of 1 . 15 ± 0.05 M© for the chemical composition of 
X = 0.55, Y = 0.30, Xcno = 0. 10, X Ne = 0.03, and Z = 0.02. The 
V and visual magnitudes deviate from our model light curve 
about 100 days after the outburst. This is due probably to the 
contribution of emission lines such as [O III] in the nebular 
ph ase. The orbital period of P or b = 0.156 days was obtained 
by Dobrotkaet al. (2006) from the modulations with a very 
small amplitude of 0.015 mag. The mass of the donor star is 
estimated to be Mi = 0.34 M© from equation (0. The cor- 
responding separation is a = 1 .4 Rq and the effective radius 
of the WD Roche lobe is R* = 0.68 R Q . We have obtained 
the epoch of e mag w e 10t to be day 49 for our 1.15 M© WD 
model as listed in Table Q] This timescale is consistent with 
the peak of the secondary maximum and also the emergence 
of the companion star from the WD envelope. 

V2362 Cyg was discovered by Nishimura at mag 10.5 on 
2006 April 2.807 UT. The nova was not detected on his films 
taken on March 28 (limiting mag 12) or on earlier patrol films 
back to 2001. Therefore, we assume here the outburst day 
of f 0B = 2453827.0 (April 1.5 UT). Optical and near infrared 
light curves are plotted in Figure [4] Our theoretical light 
curves of free-free emission again reasonably fit with the / 
observation until day ~ 500 except the secondary maximum. 
The best-fit model is the WD mass of 0.7 ±0.05 M© for the 
chemical composition of X = 0.45, Y = .18, Xcno = 0.35, 
and Z = 0.02 (see, e.g.. lMunari et al. Il2008l for observed val- 
ues). The R and y magnitudes slightly but the V and visual 
magnitudes largely deviate from our model light curve after 
the secondary maxim um. The orbital period of P or \, = 0.207 
days was obtained by Goranskij et al. (2008) from the mod- 
ulations with an amplitude of 0. 1 1 mag. If we take the donor 
mass of M2 = 0.48 M© from equation (0, the separation is 



a = 1.6 R Q , and the WD Roche lobe is R* x = 0.66 Rq- We 
have examined the epoch of e mag « e rot to be day 240 for our 
model of 0.7 Mq WD as listed in Table Q] This timescale 
is again reasonably consistent with the peak of the secondary 
maximum and also the emergence of the companion star from 
the WD envelope. 

5. CONCLUSIONS 

We have estimated the WD mass of the classical novae 
V2491 Cyg by comparing our free-free light curves with the 
optical and near infrared observations as well as by compar- 
ing our blackbody X-ray light curves with the Swift XRT data. 
The best-fit model is the 1 .3 ± 0.02 M© WD for the chemical 
composition of X = 0.20, Y = 0.48, X CN0 = 0.20, X Ne = 0.10, 
and Z = 0.02. We strongly recommend composition analysis 
of the ejecta because the enrichment of WD matter provides 
information whether the WD mass increases or not. We have 
also estimated the WD mass of V1493 Aql and V2362 Cyg 
to be M W d = 1.15 ±0.05 and 0.7 ±0.05 Mq, respectively. 
For these three novae, the epoch of secondary maximum is 
consistently explained by our magnetic activity model if the 
magnetic activity reaches maximum at e mag ~ e ro t in the WD 
envelope. We strongly recommend search for magnetic activ- 
ities for these three novae even in quiescence. 

We thank D. Takei for providing us with their machine read- 
able X-ray data of V2491 Cyg and also AAVSO and VSOLJ 
for the optical and near infrared data for V2491 Cyg, V1493 
Aql, and V2362 Cyg. We are grateful to the anonymous ref- 
eree for useful comments and to D. Takei, M. Tsujimoto, and 
S. Kitamoto for stimulating discussion on V2491 Cyg. This 
research has been supported in part by the Grant-in-Aid for 
Scientific Research (20540227) of the Japan Society for the 
Promotion of Science. 
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